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^1^, Abstract. Ionized gas and stellar kinematics have been measured along the major axes of seventeen nearby spiral 

■ galaxies of intermediate to late morphological type. We discuss the properties of each sample galaxy distinguishing 
between those characterized by regular or peculiar kinematics. In most of the observed galaxies ionized gas rotates 
more rapidly than stars and have a lower velocity dispersion, as is to be expected if the gas is confined in the 
disc and supported by rotation while the stars are mostly supported by dynamical pressure. In a few objects, 

I I gas and stars show almost the same rotational velocity and low velocity dispersion, suggesting that their motion 

^ . is dominated by rotation. Incorporating the spiral galaxies studied by Bertola et al. (1996), Corsini et al. (1999, 

QQ ' 2003) and Vega Beltran et al. (2001) we have compiled a sample of 50 SO/a-Scd galaxies, for which the major-axis 

1/^ , kinematics of the ionized gas and stars have been obtained with the same spatial (~ 1") and spectral (~ 50 

1/^ ■ km s~^) resolution, and measured with the same analysis techniques. This allowed us to address the frequency of 

' counterrotation in spiral galaxies. It turns out that less than 12% and less than 8% (at the 95% confidence level) of 

, the sample galaxies host a counterrotating gaseous and stellar disc, respectively. The comparison with SO galaxies 

JnI ■ suggests that the retrograde acquisition of small amounts of external gas gives rise to counterrotating gaseous 

I discs only in gas-poor SOs, while in gas-rich spirals the newly acquired gas is swept away by the pre-existing 

. gas. Counterrotating gaseous and stellar discs in spirals are formed only from the retrograde acquisition of large 

CL(' amounts of gas exceeding that of pre-existing gas, and subsequent star formation, respectively. 

2 ' Key words, galaxies: kinematics and dynamics - galaxies: spiral - galaxies: structure 

-4— » ■ 

• 

^ 1. Introduction ences therein) , and the presence of pressure-supported ion- 

• . ized gas in bulges (Bertola et al. 1995; Cinzano et al. 1999). 

^ ■ Studying the interplay between ionized gas and stellar All these issues will greatly benefit from a survey de- 

^ . kinematics allows us to address different topics concerning voted to the comparative measurements of ionized gas and 

■ - - ■ the dynamical structure of disc galaxies and to constrain stellar kinematics in SOs and spiral galaxies. Since these 

the processes leading to their formation and evolution, r^^e available only for a limited number of objects, in the 

These topics include the study of the mass distribution of past years we began a scientific programme aimed at de- 

luminous and dark matter (see Sofue & Rubin 2001 for living detailed velocity curves and velocity dispersion ra- 

a review), the ubiquity of supermassive black holes and dial profiles of ionized gas and stars along the major axes 

their relationship with the large-scale properties of the of disc galaxies (Bertola et al. 1995, 1996; Corsini et al. 

host galaxies (see Merritt & Ferrarese 2001 for a review), 1999^ 2003; Vega Beltran et al. 2001; Funes et al. 2002) to 

the discovery of kinematically decoupled components (see be used for mass modeUing (Corsini et al. 1999; Cinzano 

Bertola & Corsini 1999 for a review), the origin of disc et al. 1999; PignateUi et al. 2001). 

heating (Merrifield, Gerssen & Kuijken 2001 and refer- xhis paper is organized as follows. An overview of 

the properties of the sample galaxies as well as the spec- 

~ ; ~ ~~ ; ~ . „ ^ , , ., troscopic observations and data analysis are presented in 

bend offprint requests to: pizzellaiSpd. astro. it . m mi ^ ■ ■ ■ ^ ^ n i • 

u i- -J iiiiT^ bectionlzl The resulting ionized gas and stellar kinematics 

Based on observations carried out at the iiuropean ' _ 

Southern Observatory La SiUa (Chile) (ESO 56.A-0684 and are given and interpreted in Section |31 In particular, we 

57 A-0569) derived the central velocity dispersion of stars of all the 

** Tables 3 and 4 are only available in electronic form at the sample galaxies with the aim of studying the relationship 

CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) between that the disc circular velocity and bulge velocity 

or via .http://cdsweb.u-strasbg.fr/Abstract.html, dispersion (Ferrarese 2002; Baes et al. 2003) in a forthcom- 
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Fig. 1. Absolute magnitude distribution for the sample 
galaxies. A line marks M^^ — —20.05, which corresponds 
to M* for spiral galaxies as derived by Marzke et al. (1998) 
and assuming Hq — 75 km Mpc~^. The dashed re- 
gion identifies galaxies classified barred or weakly barred 
in RC3. 

ing paper (but see also Pizzella et al. 2003). The fraction 
of spiral galaxies hosting a counterrotating component is 
estimated in Section 0] by analysing the major-axis kine- 
matics of all the spiral galaxies we have observed in recent 
years. Our conclusions are discussed in Section |5| 

2. Sample selection, spectroscopic observations 
and data reduction 

2. 1 . Sample selection 

All the observed galaxies are bright {Bt < 13.5) and 
nearby {Vq < 5900 km s^^) with an intermediate-to- 
high inclination (40° < i < 85°). The Hubble morpho- 
logical type of the sample galaxies ranges from Sb to Scd 
(de Vaucouleurs et al. 1991, RC3 hereafter) with eight 
unbarred, seven weakly barred and two strongly barred 
galaxies (RC3). The galaxies have been chosen for their 
strong emission lines. An overview of the basic properties 
of the sample galaxies is given in Table ^ The distribu- 
tion of their absolute magnitudes, which brackets the M* 
value for spiral galaxies (Marzke et al. 1998, for Hq = 75 
km s~^ Mpc~^) is shown in Figure 

2.2. Spectroscopic observations 

The spectroscopic observations of our sample galaxies 
were carried out at the ESO in La Silla with the 1.52 m 
ESO telescope on 1996 January 22-27 (run 1 hereafter), 
and on 1996 August 9-14 (run 2 hereafter). 



The telescope mounted the Boiler & Chievens 
Spectrograph. The grism No. 26 with 1200 grooves mm~^ 
blazed at 5730 A was used in the first order in combi- 
nation with the 2'.'5 x 4^2 slit and the Loral CCD No. 
24 with 2048 x 2048 pixels of 15 x 15 ^m^. The wave- 
length range between 4940 A and 6940 A was covered with 
a reciprocal dispersion of 0.98 A pixel"^, which guaran- 
tees adequate oversampling of the instrumental broaden- 
ing function. Indeed, the instrumental resolution, obtained 
by measuring the width of emission lines of a compari- 
son spectrum after the wavelength calibration, was 2.59 A 
(FWHM). This corresponds to an instrumental velocity 
dispersion of ct = 1.10 A (i.e. ~ 70 and « 50 km s~^ at 
the blue and red edges of the spectra, respectively). The 
angular sampling was 0'.'81 pixel""'^. 

At the beginning of each exposure, the slit was cen- 
tred on the galaxy nucleus using the guiding TV cam- 
era and aligned along the galaxy major axis. Details of 
the slit position and spectra exposure times are given in 
Table|5| Comparison lamp exposures were obtained before 
and/or after each object integration to allow an accurate 
wavelength calibration. Quartz lamp and twilight sky fiat- 
fields were used to map pixel-to-pixel sensitivity variations 
and large scale illumination patterns. A dozen of late G 
and early K stars were observed with the same set up 
to serve as templates in measuring the stellar kinematics. 
The value of the seeing FWHM during the different ob- 
serving runs ranged between I'.'O and 2'.'0, as measured by 
fitting a two-dimensional Gaussian to the guide star. 

Table 2. Log of spectroscopic observations 



Object Run tcxp P.A. 



[name] [s] 



NGC 210 


2 


5 


X 


3600 


160 


NGC 615 


2 


7 


X 


3600 


155 


NGC 1620 


1 


7 


X 


3600 


25 


NGC 2590 


1 


4 


X 


3600 


77 


NGC 2708 


1 


5 


X 


3600 


20 


NGC 2815 


1 


5 


X 


3600 


10 


NGC 3054 


1 


5 


X 


3600 


118 


NGC 3200 


1 


4 


X 


3600 


169 


NGC 3717 


1 


3 


X 


3600 


33 


NGC 4682 


1 


3 


X 


3600 


85 


NGC 5530 


2 


3 


X 


3600 


127 


NGC 6118 


2 


5 


X 


3600 


58 


NGC 6878 


2 


6 


X 


3600 


125 


NGC 6925 


2 


8 


X 


3600 


5 


NGC 7083 


2 


6 


X 


3600 


5 


NGC 7412 


2 


6 


X 


3600 


65 


NGC 7531 


2 


5 


X 


3600 


15 



2.3. Routine data reduction 

The spectra were bias subtracted, flat-field corrected, 
cleaned for cosmic rays and wavelength calibrated using 
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Table 1. Parameters of the sample galaxies 



Object 




Type 


Bt 


PA 


i 


v@ 


D 


Scale 


R25 




-Rlast 

[stars] 


/-R25 
[gas] 


[name] 


[RSA] 


[RC3] 


[mag] 


n 


[°] 


n — li 

[km s ] 


[Mpc] 


r ///I 

[pc/ ] 


[] 


[mag] 












(4.) 






y ' I 












yio) 


NGC 210 


Sb(rs) 


SABb(s) 


11.60 


165 


49 


1650 


23.2 


113 


150 


-20.5 


0.3 


0.8 


NGC 615 


Sb(r) 


SAb(rs) 


12.47 


155 


67 


1849 


25.7 


125 


108 


-20.3 


0.7 


0.9 


NGC 1620 




SABbc(rs) 


13.08 


25 


70 


3509 


46.1 


224 


86 


-21.1 


0.7 


1.2 


NGC 2590 




SAbc(s): 


13.94 


77 


72 


4960 


63.4 


308 


67 


-21.0 


0.5 


1.0 


NGC 2708 




SABb(s)pec? 


12.80 


25 


60 


1984 


23.5 


114 


78 


-19.4 


0.6 


0.8 


NGC 2815 


Sb(s) 


SBb(r): 


12.81 


10 


72 


2535 


30.0 


145 


104 


-21.0 


0.4 


1.0 


NGC 3054 


SBbc(s) 


SABb(r) 


12.35 


118 


52 


2430 


28.5 


138 


114 


-20.3 


0.5 


0.9 


NGC 3200 


Sb(r) 


SABc(rs): 


12.83 


169 


73 


3526 


43.4 


211 


125 


-21.5 


0.5 


0.9 


NGC 3717 


Sb(s) 


SAb:sp 


12.24 


33 


81 


1748 


19.6 


95 


180 


-20.5 


0.4 


0.5 


NGC 4682 


Sc(s) 


SABcd(s) 


13.14 


83 


61 


2344 


28.8 


140 


77 


-20.0 


0.3 


0.8 


NGC 5530 


Sc(s) 


SAbc(rs) 


11.79 


127 


62 


1092 


11.7 


57 


125 


-19.4 


0.6 


0.9 


NGC 6118 


Sc(s) 


SAcd(s) 


12.42 


58 


65 


1580 


21.4 


104 


140 


-20.4 


0.6 


0.8 


NGC 6878 


Sc(r) 


SAb(s) 


13.45 


125 


40 


5821 


77.3 


375 


48 


-21.3 


0.6 


1.0 


NGC 6925 


Sbc(r) 


SAbc(s) 


12.07 


5 


75 


2783 


37.7 


183 


134 


-21.9 


0.6 


0.9 


NGC 7083 


Sb(s) 


SAbc(s) 


11.87 


5 


53 


3093 


39.7 


192 


116 


-21.5 


0.5 


1.0 


NGC 7412 


Sc(rs) 


SBc(s) 


11.88 


65 


42 


1712 


22.7 


110 


116 


-20.1 


0.4 


0.8 


NGC 7531 


Sbc(r) 


SABbc(r) 


12.04 


15 


67 


1591 


20.9 


102 


134 


-20.2 


0.4 


0.6 



NOTES - ^■^•'Morphological classification from Sandage & Tammann (1981, RSA hereafter). '''-'Morphological classification from 
RC3. Total observed blue magnitude from RC3 except for NGC 3954, NGC 4682 and NGC 5530 (LEDA). Major-axis 
position angle from RC3. '^■'inclination, derived as cos^ i = {q^ — qo)/(l ~ 1o)- The observed axial ratio q is taken from RC3 
and the intrinsic flattening qo = 0.11 has been assumed following Guthrie (1992). ''^^ Heliocentric velocity of the galaxy derived 
at the centre of symmetry of the rotation curve of the gas. AV© = 10 km s~^. Distance obtained as Vo/Ho with Ho = 75 
km s~^ Mpc~^ and Vo, the systemic velocity derived from Vq corrected for the motion of the Sun with respect to the Local 
Group, as in RSA. Radius of the 25 B mag arcsec ^ isophote derived as R25 = -D25/2 with D25 from RC3. '^^^ Absolute total 
blue magnitude corrected for inclination and extinction from RC3. '^^^ Radial extension of the stellar rotation curve in units of 
R25 from this paper. ''^^^ Radial extension of the ionized gas rotation curve in units of R25 from this paper. 



standard MIDAS"'^ routines. Cosmic rays were identified by 
comparing the counts in each pixel with the local mean 
and standard deviation (as obtained from the Poisson 
statistics of the photons knowing the gain and readout 
noise of the detector), and then corrected by interpolat- 
ing a suitable value. In the central regions (typically for 
jrj ^ 5") cosmic rays have been removed by manually 
editing the spectra. 

The instrumental resolution as a function of wave- 
length was derived as the Gaussian FWHMs measured 
for several unblended arc lamp lines distributed over the 
whole spectral range of a wavelength-calibrated compar- 
ison spcctriim. Finally, the spectra of the same galaxy 
were aligned and co-added using the centers of their stel- 
lar continua as reference. In the resulting spectra, the 
contribution from the sky was determined by interpolat- 
ing along the outermost 10"-20" at the edges of the slit, 
where galaxy light was negligible. The sky level was then 
subtracted. Since the spatial extension of the emission 
lines can exceed that of the stellar continuum, different 
ranges for sky subtraction have been considered in ion- 



^ MIDAS is developed and maintained by the European 
Southern Observatory. 



ized gas and stellar kinematic measurements. In few cases 
the galaxy emission lines extended through the whole slit 
length. The sky level was derived by scaling the sky frame 
derived for a different galaxy to match the intensity of the 
relevant night sky emission lines. 

2.4. Measuring ionized gas and stellar kinematics 

The ionized gas kinematics was measured by the simul- 
taneous Gaussian fit of the emission lines present in the 
spectra (namely [N II] A6583, Ha, and [S II] AA6716, 6731). 
The galaxy continuum was removed from the spectra, as 
was done for measuring the stellar kinematics. We fit- 
ted a Gaussian to each emission line in each row of the 
continuum-subtracted spectrum, assuming them to have 
the same line-of-sight velocity (Vg), and velocity dispersion 
(cTg). Velocities and velocity dispersions were corrected 
for heliocentric velocity and instrumental FWHM, re- 
spcictively. An additional absorption Gaussian was added 
to the fit to take into account for the presence of the 
Ha absorption line and the flux ratio of the [N II] lines 
have been fixed at 1:3. Far from the galaxy centre (for 
|r| ^ 10") we averaged adjacent spectral rows in order to 
increase the signal-to-noise ratio of the relevant emission 
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lines. We checked for a few galaxies that the error in the 
kinematic parameter determination derived by means of 
Monte Carlo simulations did not differ significantly from 
the formal errors given as output by the least-squares fit- 
ting routine. We therefore decided to assume the latter as 
error bars on the gas kinematics. For each galaxy we de- 
rive the heliocentric systemic velocity, Vq, as the velocity 
of the centre of symmetry of the ionized gas rotation curve 
(Tabled. 

The stellar kinematics were measured from the galaxy 
absorption features present in the wavelength range run- 
ning from 5050 A to 5550 A and centred on the Mg line 
triplet (AA 5164, 5173, 5184 A). We used the Fourier cor- 
relation quotient method (FCQ, Bender 1990) following 
the prescriptions of Bender, Saglia & Gerhard (1994). The 
spectra were binned along the spatial direction to obtain a 
nearly constant signal-to-noise ratio larger than 20 per res- 
olution element. The galaxy continuum was removed row 
by row by fitting a fourth to sixth order polynomial as in 
Bender et al. (1994). The Kl III star HR 5777 was adopted 
as a kinematic template. This allowed us to derive, for each 
spectrum, the line-of-sight stellar velocity (?;*), and veloc- 
ity dispersion (fi) by fitting a Gaussian to the line-of-sight 
velocity distribution (LOSVD) at each radius. Velocities 
were corrected for heliocentric velocity. We derived errors 
for the stellar kinematics from photon statistics and CCD 
read-out noise, calibrating them by Monte Carlo simula- 
tions as done by Gerhard et al. (1998). These errors do 
not take into account possible systematic effects due to 
any template mismatch. 

The ionized gas and stellar kinematics of all the sample 
galaxies are tabulated in Tables 3 and 4, respectively, and 
plotted in Figure |21 

3. The ionized gas and stellar kinematics 

3.1. Velocity curves and velocity dispersion profiles of 
ionized gas and stars 

The resulting ionized gas and stellar kinematics of all our 
sample galaxies are shown in Figure|21 For each object the 
plot is organized as it follows: 

1. In the upper panel we display the galaxy image ob- 
tained from the Digitized Sky Survey. The galaxy im- 
age has been rotated and the slit position has been 
plotted to allow a better comparison between morpho- 
logical and kinematic properties. The slit width and 
length correspond to those adopted in obtaining the 
spectra. 

2. In the middle panel we plot the velocity curves of the 
gaseous (open circles) and stellar (filled circles) com- 
ponents. The velocity scale on the left of the plot in- 
dicates the observed velocity after subtraction of the 
systemic velocity given in Tabled The velocity scale 
on the right indicates the rotation velocity after the de- 
projection for the galaxy inclination given in Tabled 
Error bars are not plotted when smaller than symbols. 
The position angle of the slit is specified. 



3. In the lower panel we plot in radial profiles the veloc- 
ity dispersion of the gaseous (open circles) and stellar 
(filled circles) components. Error bars are not plotted 
when smaller than their symbols. 

The kinematic measurements obtained for the sample 
galaxies extend typically to « 0.9 R25 for the gaseous 
and « 0.5 R25 for the stellar component. The curves are 
generally regular and symmetric about the centre. 




Fig. 2. (continue). 



3.2. Comparison with the literature 

The ionized gas and/or stellar kinematics of some of the 
sample galaxies have been measured by other authors. We 
used these data sets to assess the accuracy and reliability 
of our measurements. In most cases differences between 
different authors are due to slit centring and/or position- 
ing, different analysis techniques, or both (see Fisher 1997 
for a discussion). 

Ionized gas velocity curves have been already measured 
along the major axes of NGC 615, NGC 1620, NGC 2708, 
NGC 2815, NGC 3054, NGC 3200, NGC 4682, NGC 6118, 
NGC 7083 and NGC 7531. However, these measurements 
had previously been obtained with a lower spatial and/or 
velocity resolution with respect to ours. In Figure|3|we plot 
our kinematic measurements compared to the available 
data. The agreement is good and the discrepancies {AV ^ 
20 km s~^) are due to the different spatial sampling and 
accuracy of velocity measurements. Only in the case of 
NGC 7083 do we notice a systematic difference {AV « 
80 km s~^) on the north-east side of the rotation curve 
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Fig. 2. Ionized gas {open circles) and stellar {filled circles) kinematics measured along the major axes of the sample 
galaxies. Error bars smaller than their symbols are not plotted. 



between our data and those of Rubin et al. (1982). We 
attribute such a discrepancy to the wavelength calibration 
of Rubin et al. (1982) since it is not present when we 
consider the data of Mathewson, Ford & Buchhorn (1992). 

In Fig. ^ we compared our kinematic measurements 
for the stellar component of NGC 615, NGC 2815 and 
NGC 3200 to the data derived by Bottema (1992), 
Sil'chenko, Vlasyuk & Alvarado (2001) and Vega Bcltran 
et al. (2001), respectively. The overall agreement is good 



both for velocities and velocity dispersions, except for 
NGC 2815. Indeed, we measured on the south-east side 
of this galaxy lower rotation velocities and higher velocity 
dispersions with respect to those of Bottema (1992). 

3.3. The central velocity dispersion of stars 

A tight correlation between the central velocity disper- 
sion of the spheroid and the galaxy circular velocity has 





recently been found for a sample of elliptical and spiral 
galaxies (Ferrarese 2002; Baes et al. 2003). Since the bulge 
velocity dispersion and circular velocity are related to the 
masses of the central black hole (Ferrarese & Merritt 2000; 
Gebhardt ct al. 2000) and dark matter halo (e.g. Bullock 
et al. 2001), respectively, the correlation between velocity 
dispersion and circular velocity is equivalent to unity be- 
tween the mass of the black hole and the mass of the dark 
halo (Ferrarese 2002). 



With the aim of exploring this correlation in a forth- 
coming paper (Pizzella et al. 2004. in preparation, but 
see also Pizzella et al. 2003), we tested how the central 
stellar velocity dispersion cr*(0) depends on the radius of 
the aperture in which it has been measured. Indeed when 
comparing the values of o-*(0) to those derived by other 
authors it is necessary to take this effect into account. 
Ferrarese (2002) took as reference aperture that corre- 
sponding to Re/ 8 where is the effective radius of the 
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Fig. 2. (continue). 



bulge, and adopted the aperture correction by Jorgensen, 
Franx & Kjaergaard (1995). 

To test the reliability this correction when applied to 
spiral galaxies, we measured the value of cr*(0) for each 
sample galaxy in different radial bins obtained by averag- 
ing the central 1,3, . . . , 31 rows of the spectrum (corre- 
sponding to X = O'.'Sl, 2'.'43, . . . , and 25'.'1, respectively). 
In Figure 13 the resulting profiles of (7^,(0) as a function 
of the radius of the circularized aperture are compared 
to that obtained for a galaxy with cr^(0)= 150 km s~^ 



and i?c — S'.'O by applying the aperture correction by 
Jorgensen et al. (1995). These are typical values for the 
central velocity dispersions and the bulge effective radii of 
our sample galaxies (Baggett, Baggett & Anderson 1998). 

It turns out that in most of cases the observed (7^(0) 
profile is not centrally peaked and is nearly constant or 
even rising for large-aperture radii. This suggests that 
the empirical correction by Jorgensen et al. (1995), which 
was originally proposed for elliptical galaxies with cen- 
trally peaked velocity dispersions but negligible rotation 





Fig. 2. (continue). 



velocities, is no longer valid for our late-type spiral galax- 
ies. In this case the slit aperture is positioned along the 
galaxy major axis; therefore, the contribution of the ro- 
tation velocity (which remains spatially unresolved after 
the binning process) to the measured velocity dispersion 
increases with radius. The contribution of the rotation ve- 
locity makes the actual orientation of the slit is a critical 
parameter when the aperture correction is applied. 

We therefore decided to not apply any aperture correc- 
tion to the measured value of (7^,(0). However, as shown in 



Figure |31 there is no strong dependence of (T*(0) on aper- 
ture for all the sample galaxies, except for NGC 4682. This 
peculiar behaviour is due to the low signal-to-noise ratio 
(« 10) for the innermost two radial bins. 
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Fig. 3. The ionized gas velocities derived in this study for NGC 615, NGC 1620, NGC 2708, NGC 2815 and NGC 3054 

compared with those obtained by other authors: MFB92 ~ Mathewson, Ford & Buchhorn (1992); R+82 ~ Rubin et 
ah (1982); SVAOl = Sil'chenko, Vlasyuk & Alvarado (2001); C97 = Courteau (1997); S+98 = Sofue et al. (1998). 



3.4. Interplay between kinematics of ionized gas and 
stars 

Building a complete dynamical model to address the ques- 
tion of the mass distribution of the sample galaxies is 
beyond the scope of this paper. However, it is possible 
to derive some hints about their structure directly from 
the analysis of the interplay between the kinematics of 
their gas and stars. At each radius, 14 (= \v* — Vq\) and 
Vg (= \vg — Vq\) are the observed rotational velocities of 
the stars and the ionized gas, respectively. Following Vega 
Beltran et al. (2001), we distinguished two classes in the 
sample galaxies, according to their kinematics and assum- 
ing that the gaseous and stellar discs are coplanar: 

(i) Galaxies in which ionized gas rotates faster than stars 
and has a lower velocity dispersion than the stars (i.e., 
Vg > K and (Tg < 13): NGC 210, NGC 615, NGC 1620, 
NGC 2590, NGC 2708 (ahhough it has asymmetric rota- 
tion curves for both gas and stars), NGC 2815, NGC 3200, 
NGC 3717, NGC 6118, NGC 6878, NGC 6895 and NGC 



7083. The different kinematic behaviour of the gaseous 
and stellar components can be explained by a model where 
the gas is confined in the disc and supported by rotation, 
while the stars are mostly supported by dynamical pres- 
sure. This is consistent with the application of the asym- 
metric drift (e.g. Binney & Tremaine 1987). 

(ii) Galaxies for which Vg ~ and ^ fi: NGC 4682, 
NGC 5530 and NGC 7412. The observed difference be- 
tween the velocity dispersion of the two components is 
Act < 50 km s~^. The motions of both ionized gas and 
stars are dominated by rotation. 

The weakly barred spirals NGC 3054 and NGC 7531 
show peculiar kinematics that are unexpected if the gas is 
moving in circular orbits in a disc coplanar to the stellar 
one. In fact, in NGC 3054 Vg < T4 and (Xg < fi for \r\ ^ 
10", and Vg > Vi, and o-g ~ 6 at larger radii, while NGC 
7531 has intermediate properties between the two classes. 
It is characterized by ~ K and ag < S> for \r\ ^ 10", 
and Vg > V^ and cxg ~ fi for |r| > 10". These cases cannot 
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Fig. 3. (continue). The ionized gas velocities derived in this study for NGC 3200, NGC 4682, NGC 6118, NGC 7083 
and NGC 7531 compared with those obtained by other authors: V+01 = Vega Beltran et al. (2001).; MFB92 = 
Mathewson, Ford & Buchhorn (1992); R+82 = Rubin et al. (1982); RTF80 = Rubin, Thonnard k Ford (1980); B87 
= Buta (1987); M84 = Meyssonnier (1984). 



be explained without invoking a warp in the gaseous disc 
or non- circular gas motions. 

4. The frequency of counterrotation in spiral 
galaxies 

4.1. Detectability of counterrotation 

The presence of counterrotation in disc galaxies is gener- 
ally interpreted as the signature of the external origin of 
the gaseous component (Bertola, Buson & Zeilingcr 1992). 
Since we are going to discuss the frequency of counterro- 
tation in disc galaxies, it is necessary to establish our ca- 
pability of detecting a counterrotating component when it 
is present. 

In order to estimate the frequency of counterrotating 
components we considered the sample of spiral galaxies 
obtained from Bertola et al. (1996; one Sa galaxy), Corsini 
et al. (1999; six Sa galaxies), Vega-Beltran et al. (2001; 
sixteen Sa-Scd galaxies. NGC 224 and NGC 3031 have 



been excluded from this analysis since the available kine- 
matic data do not extend to the disc region), Corsini et 
al. (2003; ten SO/a-Sa galaxies), and this paper (seventeen 
Sb-Sc galaxies). The sample consists of 50 spiral galaxies 
(0 < T < 6), with apparent magnitude 10.1 < Bt < 13.9, 
inclination i > 35° and systemic velocity Vq < 6000 
km s~^. Their major-axis ionized gas and stellar kine- 
matics have been obtained with the same spatial (FWHM 
w 1") and spectral (dinstr ~ 50 km s^^) resolution, and 
measured with the same analysis technique. These 50 spi- 
ral galaxies represent our qualified sample, for which we 
are confident that we are able to detect any counterro- 
tating or kinematically decoupled gaseous and/or stellar 
component. 

Revealing the presence of ionized gas that is counter- 
rotating with respect to the stellar component is straight- 
forward. Indeed in most a cases visual inspection of the 
shape and orientation of spectral emission and absorption 
lines enables us to detect when the gas and stars are rotat- 
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Fig. 4. The stellar kinematics derived in this study for NGC 615, NGC 2815 and NGC 3200 compared with those 
obtained by other authors: B92 = Bottema (1992); SVAOl = Sil'chenko, Vlasyuk & Alvarado (2001); V+01 = Vega 
Beltran et al. (2001). 



ing opposite directions (e.g. NGC 4450; Rubin, Graham 
& Kenney 1992). The techniques we adopted to derive the 
gaseous and stellar kinematics allowed us to measure dif- 
ferences of a few km s"-'^ in the rotation velocities of the 
two components. We are therefore confident that we are 
able to detect all the counterrotating gaseous components 
hosted in our sample galaxies. 

Disentangling a counterrotating stellar component is a 
more difficult task. In fact, the kinematics of two coun- 
terrotating stellar populations is measured from the same 
absorption features and depends on the dynamical status 
of the components, as well as on the spectral resolution 
of the available data. We estimated the upper limit on 
the fraction of counterrotating stars we can detect assum- 
ing that the counterrotating components are discs with 
the same scale length and different luminosity as follows. 
We assumed the stars in the disc of NGC 615 as the pro- 
grade component of the model. We selected the spectrum 
of NGC 615 since the absorption lines in the wavelength 



and radial ranges where we are going to measure the stel- 
lar kinematics are nearly symmetric with respect to the 
centre, have a large radial extension and are characterized 
by a high signal to noise (S/N ^ 30 per A for \r\ > 20", 
and S/N > 50 per A for |r| < 20"). We obtained the ret- 
rograde component of model by flipping the spectrum of 
NGC 615 about the centre of the slit. This corresponds 
to inverting all the line-of-sight velocities of the prograde 
component with respect to the systemic velocity of the 
galaxy. We built a series of synthetic spectra with a differ- 
ent fraction of counterrotating stars (/cr = 5%, 10%, . . . , 
50%) by means of a linear superposition of the prograde 
and retrograde models. The retrograde model has been 
divided by its own continuum and then multiplied by the 
continuum of the prograde model to remove asymmetries 
in the continuum shape on the two sides of the galaxy disc. 
We measured the stellar kinematics in the disc-dominated 
region and derived the radial profile of LOSVD of the syn- 
thetic spectra by applying the FCQ method as described 
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Fig. 5. Percentage variation in the central velocity disper- 
sion, (To , of the sample galaxies as a function of the radial 
bin, X (or of the circularized aperture, rap) in which it has 
been measured. The dotted line is the percentage varia- 
tion obtained for an elliptical galaxy with <7h.(0) — 150 
km s~^ and i?c = 8" by applying the empirical correction 
by Jorgensen et al. (1995) . 



in Sect ion [T^ We recognized the signature of the presence 
of two counterrotating components, namely the splitting 
of the LOSVD into two peaks, for /cr > 15%. We assume 
this value to be our limiting sensitivity for the detection 
of counterrotating stars in the galaxy disc. We did not 
consider the bulge-dominated region, where disentangling 
the fraction of counterrotating stars requires further con- 
straints on the dynamical status of the bulge. We conclude 
that no more than 15% of the stars in the disc of our sam- 
ple galaxies reside in retrograde orbits. 

A similar analysis has been done by Kuijken, Fisher & 
Merrifield (1996) who found an upper limit of 5% on coun- 
terrotating stars. We attribute this difference to the fact 
that they created the prograde and retrograde models by 
convolving the spectrum of a kinematic template with the 
appropriate LOSVD while we adopted a real (and there- 
fore more noisy) spectrum. It has to be noted that this is 
a limit on the relative surface brightness of the prograde 
and retrograde component. If the retrograde component 
is more radially concentrated than the prograde one (e.g. 



NGC 3593; Bertola et al. 1996), it could be detected even 
if its luminosity is less than 15% of the galaxy luminosity. 

4.2. Fraction of disc galaxies hosting a counterrotating 
gas component 

We found that ionized gas is counterrotating with respect 
to stars in 2 out 50 galaxies. They are the Sa NGC 3593 
(Bertola et al. 1996) and SO/Sa NGC 7377 (Gorsini et al. 
2003), which corresponds to a frequency of 4%. Applying 
Poisson distribution statistics, this means that less than 
12% (at the 95% confidence level) of the sample galax- 
ies have a counterrotating gaseous component. Only a few 
other spirals are known to host a counterrotating gaseous 
disc (NGC 3626, Giri, Bettoni & Galletta 1996; Haynes 
et al. 2000; NGC 4138, Jore, Broeils & Haynes 1996; 
Haynes et al. 2000; NGC 7217, Merrifield & Kuijken 1994) 
or a kinematically decoupled gas component (NGC 5854, 
Haynes et al. 2000; and the spirals hosting an inner polar 
disc listed by Gorsini et al. 2003). 

For comparison, a similar analysis can be done for the 
lenticular galaxies. Combining the data of Bertola et al. 
(1992, 1995), Kuijken et al. (1996), and Kannappan & 
Fabricant (2001) with those for the SO NGC 980 of Vega 
Bcltran et al. (2001), we obtained a sample of 53 SO galax- 
ies for which the kinematics of the ionized gas and stel- 
lar components have been measured along (at least) the 
major axis. A counterrotating or kinematically decoupled 
gaseous component has been found found in seventeen 
of these galaxies. We also took into account the galaxies 
whose ionized gas angular momentum is misaligned with 
respect to that of the stellar component since we aimed 
at deriving the fraction of SO galaxies hosting a gaseous 
component of external origin. This fraction corresponds to 
32%^J^ (at the 95% confidence level) from Poisson distri- 
bution statistics and is consistent with previous results by 
Bertola et al. (1992), Kuijken et al. (1996) and Kannappan 
& Fabricant (2001). 

4.3. Fraction of disc galaxies hosting a counterrotating 
stellar component 

In our qualified sample of 50 spiral galaxies we found 
the presence of a counterrotating stellar disc in the Sa 
NGC 3593 (Bertola et al. 1996). Applying Poisson dis- 
tribution statistics this means that less that 8% (at the 
95% confidence level) of the sample galaxies have a coun- 
terrotating stellar component. Only two other spirals are 
known to host two counterrotating stellar discs, one of 
these is corotating with the gaseous disc as observed in 
NGC 3593. They are the Sa NGC 4138 (Jore et al. 1996; 
Haynes et al. 2000) and the Sab NGC 7217 (Merrifield 
& Kuijken 1994). A kinematically decoupled stellar com- 
ponent orthogonally rotating with respect to the galaxy 
disc, has been observed in the bulge of the early-type spi- 
rals NGC 4698 (Gorsini et al. 1999; Bertola et al. 1999) 
and NGC 4672 (Sarzi et al. 2000). 
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Kuijken et al. (1996) estimated that less than 10% (at 
the 95% confidence level) of their SO galaxies host a sig- 
nificant fraction of countcrrotating stars. This is actually 
the case for NGC 4550 (Rubin et al. 1992; Rix et al. 1992) 
where half of the disc stars are moving in retrograde or- 
bits. 

5. Summary and conclusions 

We have measured the ionized gas and stellar kinematics 
along the major axes of seventeen intermediate- to late- 
type spiral galaxies. The rotation curves and velocity dis- 
persion profiles of ionized gas and star typically extend 
out to ~ 0.9ii25 and « 0.5ii25, respectively. 

In most of cases the different kinematic behaviour of 
ionized gas and stars can be easily explained if gas is con- 
fined in the disc and supported by rotation while the stars 
belong mostly to the bulge and are supported by dynam- 
ical pressure. However, kinematic peculiarities have been 
observed at least in NGC 3054 and NGC 7531. 

In addition, we discussed the frequency of counterro- 
tation in disc galaxies. We considered our seventeen spiral 
galaxies with those studied in previous papers (Bertola et 
al. 1996; Corsini et al. 1999, 2003; Vega Beltran et al. 2001) 
to build a qualified sample of 50 bright and nearby spirals, 
ranging from SO/a to Scd, for which the ionized gas and 
stellar kinematics have been measured along the major 
axis with the same analytical technique. We found that 
less than 12% and less than 8% (at the 95% confidence 
level) of these galaxies host a countcrrotating gaseous and 
stellar disc, respectively. For comparison, wc found that 
^ 30% of SOs host a countcrrotating gaseous disc, and 
Kuijken et al. (1996) estimated that less than 10% (at 
95% confidence level) host a significant fraction of coun- 
tcrrotating stars. 

To interpret the observed frc^qucncics of the gaseous 
and stellar countcrrotating components of disc galaxies, 
wc suggest a scenario in which SO and spiral galaxies are 
subject to external gas acquisiticm with equal probability. 

The retrograde acquisition of small amounts of exter- 
nal gas gives rise to countcrrotating gaseous discs only 
in SO galaxies, since in spiral galaxies the acquired gas 
is swept away by the pre-existing gas. The formation of 
countcrrotating gaseous discs is favoured in SO galaxies 
since they gas-poor systems, while spiral discs host large 
amounts of gas (Roberts & Haynes 1994; Bettoni, Galletta 
& Garcia-Burillo 2003), which is corotating with the stel- 
lar component. When they acquire external gas in retro- 
grade orbits, the gas clouds of the new retrograde and pre- 
existing prograde components collide, lose their centrifu- 
gal support, and accrete toward the galaxy centre. A coun- 
tcrrotating gaseous disc will be observed only if the mass 
of the newly supplied gas exceeds that of the pre-existing 
one (Lovelace & Chou 1996; Thakar & Ryden 1998). A 
countcrrotating stellar disc is the end-result of star for- 
mation in the countcrrotating gas disc. For this reason we 
observe a larger fraction of countcrrotating gaseous disks 
in SOs than in spirals. This also explains why the mass 



of countcrrotating gas in most SO galaxies is small com- 
pared to that of the stellar counterrotating components 
(Kuijken et al. 1996). The fraction of SOs with a kinemat- 
ically decoupled gas disc is consistent with the 50% that 
we expect if all the gas in SOs is of external origin (Bertola 
et al. 1992). 

Counterrotating gaseous and stellar discs in spirals are 
both the results of retrograde acquisition of large amounts 
of gas, and they arc observed with the same frequency. In 
this framework stellar counterrotation is the end result of 
star formation in the a counterrotating gaseous disc. 
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